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Signals of spinodal hadronization: strangeness trapping 
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If the deconfinement phase transformation of strongly interacting matter is of first-order and the 
expanding chromodynamic matter created in a high-energy nuclear collision enters the corresponding 
region of phase coexistence, a spinodal phase separation might occur. The matter would then 
condense into a number of separate blobs, each having a particular net strangeness that would 
remain approximately conserved during the further evolution. We investigate the effect that such 
strangeness trapping may have on strangeness-related hadronic observables. The kaon multiplicity 
fluctuations are significantly enhanced and thus provide a possible tool for probing the nature of 
the phase transition experimentally. 

PACS numbers: 25.75.-q, 25.75.Nq, 25.75.Gz 



I. INTRODUCTION 

One of the major goals of high-energy heavy-ion re- 
search is to explore the equation of state of strongly in- 
teracting matter, particularly its phase structure Q. De- 
pending on the beam energy, various regions of tempera- 
ture and baryon density can be explored. Thus systems 
with a very small net baryon density are formed at RHIC 
(y/s ~ 200 AGcV), while it is expected that the cre- 
ation of the highest possible baryon densities occurs at 
more moderate beam energies (y/s ~ 10 AGcV), such as 
those becoming available at the planned FAIR |3| . 

Our understanding of the QCD phase diagram is best 
developed at vanishing chemical potential, fig = 0, where 
lattice QCD calculations are most easily carried out. The 
most recent results indicate that the transformation from 
a low-entropy hadron resonance gas to a high-entropy 
quark-gluon plasma occurs smoothly at the temperature 
is raised, with no real phase transition being present 
On the other hand, at zero temperature most models pre- 
dict the occurrence of a first-order phase transition when 
the density is raised [ J5\, though no firm results are yet 
available for the corresponding value of the chemical po- 
tential, jUo. However, if the T = transformation is in 
fact of first order, one would expect the phase bound- 
ary to extend into the region of finite temperature and 
terminate at a certain critical endpoint, (/i c ,T c ) In- 
deed, recent lattice QCD results Q suggest the presence 
of such a first-order phase transition line and an associ- 
ated critical end-point, though its precise location is not 
well determined. 

It is therefore important to consider how this key is- 
sue could be elucidated on the basis of experimental 
data. Generally, one might expect that if the expanding 
matter created in a high-energy nuclear collision crosses 
a first-order phase-transition line then the associated 
non-monotonic behavior of the thermodynamic potential 
might have observational consequences. 

From this perspective, the enhancements of the K/tt 
ratio reported for beam energies of 20 — 30 AGeV at the 
SPS appears intriguing. Since these data present the 
only non-monotonic behavior seen so far in high-energy 



heavy-ion collisions, it appears appropriate and timely 
to study the consequences of a possible first-order phase 
transition on the production of kaons or, more generally, 
strange hadrons. 

A universal feature of first-order phase transitions is 
the occurrence of spinodal decomposition, which results 
from the convex anomaly in the associated thermody- 
namic potential jlCj . This phenomenon occurs when bulk 
matter, by a sudden expansion or cooling, is brought 
into the convex region of phase coexistence. Since such a 
configuration is thermodynamically unfavorable and me- 
chanically unstable, the uniform system seeks to reorga- 
nize itself into spatially separate single-phase domains. 
Moreover, since this spinodal phase separation develops 
by means of the most unstable collective modes, the re- 
sulting domain pattern tends to have a scale character- 
istic of those modes. This general phenomenon, which is 
known in many areas of physics and has found a variety 
of technological applications, appears to be an important 
mechanism behind the multifragmentati on p henomenon 
in medium-energy nuclear collisions 0, llOj . where the 
relevant first-order phase-transition is between the nu- 
clear liquid and a gas of nucleons and light fragments. 
Thus, if a first-order phase transition is encountered dur- 
ing the expansion stage of a high-energy nuclear collision, 
one might expect that such a spinodal separation might 
occur. While the resulting enhancement of baryon fluc- 
tuations was studied by Bower and Gavin and the 
prospects for observing such a process via iV-body kine- 
matic correlations was discussed by Randrup |l2j . the 
present study explores the consequences for the produc- 
tion of strange hadrons. 

In order for a spinodal decomposition to occur, sev- 
eral conditions must be met. First of all, of course, the 
equation of state must have a first-order phase transi- 
tion. Though expected, the existence of a first-order 
phase transition is not yet well established theoretically 
and it may ultimately have to be determined experimen- 
tally by analyzes of the kind considered here. Second, 
the dynamical trajectory of the bulk matter formed early 
on must pass through the spinodal region of phase coex- 
istence. While some calculations suggest that this may 
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happen at FAIR [l3lll3 . Il5j . this question needs to be in- 
vestigated more thoroughly. Third, even if the the above 
conditions are met, the dynamical conditions of the colli- 
sion must be carefully tuned to ensure, on the one hand, 
that the bulk of the system is brought into the spinodal 
region sufficiently quickly to achieve a quench, yet, on the 
other hand, the overall expansion should be slowed down 
to a degree that will allow the dominant instabilities to 
grow sufficiently to cause the bulk to break up. While the 
conditions for achieving this delicate balance are hard to 
ascertain theoretically, they may be found by a system- 
atic variation of the beam energy. These open questions 
notwithstanding, we shall here assume that the matter 
created in a heavy-ion collision somehow breaks into a 
number of subsystems, blobs, which subsequently expand 
and hadronizc independently and we then investigate the 
consequences for the production of strange hadrons. In 
particular, we wish to ascertain whether such a breakup 
could lead to an enhancement of the K/tt ratio and its 
fluctuations. 

In such a scenario, if the breakup is sufficiently rapid, 
then whatever net strangeness happens to reside within 
the region of the plasma that forms a given blob will 
effectively become trapped and, consequently, the result- 
ing hadronization of the blob will be subject to a corre- 
sponding constraint on the net strangeness. As we shall 
demonstrate, this type of canonical constraint will en- 
hance the multiplicity of strangeness-carrying hadrons, 
as compared to the conventional (grand-canonical) sce- 
nario where global chemical equilibrium is maintained 
through the hadronic freeze-out Oil. (This occurrence 
of an enhancement is qualitatively easy to understand, 
since the presence of a finite amount of strangeness in 
the hadronizing blob enforces the production of a cor- 
responding minimum number of strange hadrons.) The 
fluctuations in the multiplicity of strange hadrons, such 
as kaons, are enhanced even more, thus offering a pos- 
sible means for the experimental exploration of the phe- 
nomenon. 

The remainder of this paper is organized as follows: 
First we introduce a suitable idealized model framework 
and develop the necessary formal tools for the required 
canonical calculations, with some of the formal manipula- 
tions being relegated to appendices. The key features are 
then brought out in a schematic scenario containing only 
charged kaons. Subsequently, we present instructive nu- 
merical results and also make a quantitative assessment 
of the importance of global strangeness conservation. We 
finally give a concluding discussion. 



II. CALCULATION AL FRAMEWORK 

In order to establish a framework for investigating the 
effect of the strangeness trapping mechanism, we adopt 
the following schematic scenario: We start by consider- 
ing the expanding system when it is still in the plasma 
phase. At this stage the system is spatially uniform and 



the strange quarks and antiquarks can be considered as 
being randomly distributed throughout the system, irre- 
spective of what the net baryon density happens to be. 
We imagine that the bulk of the expanding and cooling 
plasma enters the region of phase coexistence and that 
the associated spinodal instability will cause it to break 
up into separate subsystems, blobs, which are assumed 
to all have the same size, as they would tend to have 
in a spinodal breakup. Each of these blobs now pro- 
ceed to expand and hadronize while maintaining its net 
strangeness. The resulting assembly of hadrons is deter- 
mined at freeze-out by a sampling of the statistical phase 
space, subject to the appropriate canonical strangeness 
constraint. 

In order to assess the effect of the strangeness trap- 
ping, it is useful to compare the results against the stan- 
dard scenario, in which the system is assumed to evolve 
to freeze-out while remaining macroscopically uniform 
and maintaining global statistical equilibrium. Focusing 
on a particular subvolume Vu, we describe the resulting 
hadron gas in the classical grand-canonical approxima- 
tion. The abundance of a particular hadron specie k is 
then 

fl 9kVhTl k ( ?^\ e ^ B B k +n Q Q k+ » s S k )/T h (1) 
2tt 2 h 3 c 3 V T h 1 ' V ; 

where K2(x) = x^Kiix) is regular at x — and a hadron 
of the specie k has baryon number B k , electric charge 
Qk, and strangeness S k . The average values of baryon 
number, charge, and strangeness in the volume Vh then 
readily follow, 

B = ^ B k n k , Q = 2J Qkfik , S = ^ S k n k • (2) 

The values of the freeze-out temperature Th and the three 
chemical potentials fig, /iq, and fig w iH be determined 
by fits to the experimental yield ratios (see later). In this 
treatment, the individual hadron species are statistically 
independent and the associated multiplicities have Pois- 
son distributions, so the multiplicity variances are equal 
to the mean values, a\ = n k . 

We now return to the particular spinodal scenario de- 
scribed above, where we assume that the plasma has bro- 
ken up into separate blobs. We first consider the distribu- 
tion of strangeness within the blobs and then treat their 
subsequent hadronic freeze-out. 

If a given plasma blob is only a small part of the total 
system, its statistical properties may be treated in the 
grand-canonical approximation. The various quark (and 
gluon) species arc then independent. Furthermore, since 
there is no bias on the overall strangeness, the s and s 
quarks have identical partition functions, Z s = Z Sl where 

\nZ s = + (3) 

3 V q Tl f ms \ 
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Here g q = 6 is the quark spin-color degeneracy, T q is the 
plasma temperature, and V q is the volume of the par- 
ticular blob considered at the time when its strangeness 
is frozen in. The energy e s is given by t 2 s = p 2 s + m s , 
where we use the mass m s = 150 MeV. Whereas the first 
expression is the exact fermionic form, the last relation 
emerges in the classical limit which we shall adopt here 
for simplicity (see Appendix |A"|) . Then the s and s mul- 
tiplicities, v and u, have Poisson distributions character- 
ized by the mean value £ s = In Z s . The total strangeness 
content in the blob is then So = v — v. Since the different 
quark flavors are distributed independently, the resulting 
probability for ending up with a given blob strangeness 
Sq is independent of the prevailing baryon and charge 
contents and can be expressed as a modified Bessel func- 
tion, 



P(S ) 



u\TA 



v-v,s = Is {2(s)e 



(5) 



We note that the corresponding ensemble average value 
of Sq vanishes, -<5o^= 0, while its ensemble variance is 

To avoid spurious correlations, we take account of the 
fact that the presence of a certain net strangeness in a 
given blob introduces a bias on its baryon number and 
charge, relative to the overall grand-canonical averages B 
and Q of the grand-canonical reference scenario, Eq. J5J). 
Indeed, each particular value of Sq determines a canonical 
subensemble of blobs that have modified distributions of 
baryon number and charge. It is elementary to show that 
these are shifted by amounts proportional to So, so the 
corresponding conditional expectation values become 



(B)s a — B - \Sq , (Q)s a — Q + kSo 



(6) 



It also follows that the ensemble correlations of B and Q 
with S are given by 



obs = <BS y = 
ctqs = <QS>- = 



So 



3 

-3°S 



(7) 
(8) 



Each isolated blob is assumed to expand further while 
hadronizing, until the freezeout volume Vh = "xVq nas 
been reached. Its temperature is then Tk .< T q . Rough 
approximations to the equation of state |l7j and the de- 
mand of energy conservation suggest that the expansion 
factor is x ~ 3, to within a factor of two or so. The 
blob has now transformed itself into a hadron resonance 
gas which we describe as a canonical ensemble charac- 
terized by the strangeness So of the precursor plasma 
blob. The baryon number and charge are treated grand- 
canonically and the demand that the expected baryon 
and charge contents match the above conditional values 
(JSJ) then determines the associated biased chemical po- 
tentials n' B and fj,Q , where the primes indicate that these 
pertain to the biased (canonical) ensemble characterized 
by the particular value of Sq. Though required for formal 



consistency and included in the calculations, this refine- 
ment is not quantitatively important. 

For the description of the hadron gas, we include 124 
hadronic species {fc}, from the 7r°(135) and up though 
the Sl — (1672). Each specie is characterized by its one- 
particle partition function, 

( k = g k J d Tk ^ P * c -(ek-n'nB k -n'nQ k )/T h ^ 

27T 2 h 3 c 3 T h 

The non-strange hadrons, which are not affected directly 
by the canonical strangeness constraint, have grand- 
canonical distributions governed by the biased chemical 
potentials [i' b (Sq) and ^!q(Sq). The total partition func- 



tion is therefore of the form Z = Z, 



_ {5=0} 7 {S#0} 



We 



describe briefly below how we obtain the canonical parti- 
tion function for the strange hadrons, Zg^°\ and refer 
to Appendix iBl for more details. 

For this task, we organize the hadron species accord- 
ing to their strangeness (an alternate approach was em- 
ployed in Refs. [S E3). For each strangeness class 
S = ±1,±2,±3, we introduce the effective one-particle 
partition function = J2kCk$s k ,s, which accounts 
for all hadron species having the particular strangeness 
S. The corresponding generic multiplicities are denoted 
by Ns, the multiplicity of hadrons having the specified 
strangeness S (e.g. N + i = Hro + n K + + + • • •)■ ^ 
should be noted that (s is not simply related to £_s when 
the system has a net baryon density. In terms of these 
quantities we then have 



Asm 
-■So 



n 



sN s 

E ss 

N s >0 ° 



5(J2n s S-Sq). (11) 



It is convenient to obtain this total partition function 
recursively by adding two conjugate strangeness classes 
{±S} at a time, 



r {S=±l,±2} 



E 

S 2 =0,±2,... 



7 {S=±1,±2,±3} 
■So 



Z {±1} z {±2} 
Sq — S2 S2 ' 



{±1,±2} -{±3} 



^ / ^So — S3 
S 3 =0,±3,... 



(12) 



(13) 



where the canonical partition function for a single pair 
of conjugate strangeness classes is given by an expression 
analogous to Eq. (J3J , 



r{±S} 



E 

N + N- 



7V+! AU (Ar+ - JV - )5 ' So 



C+s 



So/ 2 



(14) 



The corresponding correlated conditional distribution of 
the generic multiplicities {Ag} is then determined. Fur- 
thermore, recursive expressions can readily be derived for 
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the associated multiplicity moments, such as (Ns)s an( i 
(NgNs'} So (see Appendix[HJ). Once we know the number 
of hadrons with a given strangeness S, N$, we can obtain 
the multiplicities {n K } of the individual strange species 
from the corresponding binomial distributions, 



P*{n K ) 



W Q" (Cs - C) 



N s —n K 



N s 



J (N s -n K )\ 



(15) 



By proceeding as described above, it is possible to treat 
the entire ensemble of possible blobs and by numerical 
simulation generate a sample of "events" consisting of the 
resulting primordial hadrons. (The importance of subse- 
quent decays is discussed later.) These "final states" can 
then be analyzed as idealized experimental data. This 
can be done both in the spinodal scenario where individ- 
ual blobs are treated canonically as well as in the grand- 
canonical reference scenario. For instructive purposes, 
we also consider a restricted canonical scenario in which 
the strangeness of a blob is always required to vanish, 
So = 0. 

The inclusion of the bias effect (see Eq. (JB)), as ex- 
pressed through \J B and (j,'q in (JTDJ , ensures that the 
ensemble correlation of baryon number and charge with 
strangeness remains the same as it were in the plasma, 
even though our treatment of the hadron production con- 
serves B and Q only on the average for each So, 



a BS = <BS> <By<S> 

= Y,(B-\So)S P{So) = 

So 

oqs = <QS> <Q>< S >- 

= £(G + %S )S P(S ) = 

So 



J2(B)s S P(So) 

So 



3 a So > 



(16) 



J2(Q)s a SoP(So 

So 



~3 U S 



(17) 



III. SCHEMATIC SCENARIO 

Before discussing the results of such numerical sim- 
ulations, it is instructive to first illustrate the effect of 
strangeness trapping in a schematic scenario. For this 
purpose, assume that plasma blobs of strangeness So 
hadronize into charged kaons only. We furthermore ig- 
nore the chemical potentials (fis is ineffective for mesons 
and fj,Q is anyway rather small). The resulting multi- 
plicity distributions are then given in terms of the (com- 
mon) one-particle partition function, 



Ck = 9K j 



d 3 rd 3 p 



-e K /T _ 



9K VhTl r> ,m K 

w ~¥J K2{ ^ ] ■ (18) 



In particular, the average kaon multiplicities resulting 
from blobs of strangeness Sq can be expressed as an 
asymptotic expansion in 1/Ck, 



(n K ±)< 



Ck ± hSo 



4Sp 2 - 1 
32Ca- 



0( 



1 

SK 



(19) 



Ignoring the relatively small effects of quantum statistics 
on the quark multiplicity distribution, we found above 
that the ensemble of values So = v — v is characterized 
by -< So y= and a s = 2( s , The ensemble average K ± 
multiplicities are then 



< n K ± y = Ck — j 



1 



(20) 



Furthermore, the corresponding expressions for the en- 
semble multiplicity (co)variances are 



'K+ 



°K~ 



ICJC + Hs + 0{ 



1 



K 



a K+K 



- = hCK-hc + o{—) . 



K 



(21) 
(22) 



It should be noted that these expressions imply that 
a 2 (riK+ —n>K- ) = 2£s = cr| n as required by the fact that 
So = n K + - rii 



■- in each blob. 



In the special case where £ s happens to equal C,k we re- 
cover the grand-canonical result, < n K ± >-= a'j (± = Ck 
and aK+K- = 0, reflecting the fact that the grand- 
canonical ensemble can be built from an ensemble of 
canonical subcnscmblcs. But generally the above results 
differ from those of the grand-canonical scenario. In par- 
ticular, when Ck < Cs = (1 + A)Ck> as tends to be the 
case due to the large degeneracy of the quarks, the kaon 
production is enhanced by the strangeness trapping, 



-< n K ± >~ 

_2 



Ck + \A, 
Ck + |ACic 
-h&CK ■ 



(23) 
(24) 
(25) 



The negative value of the covariance reflects the fact pos- 
itive (negative) values of So favors positive (negative) 
kaons, so an excess of positive kaons is likely to be ac- 
companied by a deficit of negative kaons, and vice versa. 
We note that while the relative increase of the average 
multiplicities is of the order of A/Ck, the relative effect 
on the fluctuations is of the order of A, thus being en- 
hanced by a factor of Ck- This basic feature suggests 
that the fluctuations are preferable to the averages for 
probing the strangeness trapping phenomenon. 



IV. RESULTS AND DISCUSSION 

We now discuss results obtained by numerical sam- 
pling of the canonical partition function described above. 
In our calculations, which serve to merely illustrate the 
effect of the strangeness trapping mechanism, we use 
the freeze-out hadron populations directly in the anal- 
ysis and make no attempt to include subsequent elec- 
troweak decays. Obviously, this complication should be 
addressed before a quantitative confrontation with data 
can be made. 
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FIG. 1: The employed freeze-out values of the temperature 
Th and the chemical potentials fiQ and us as functions of 
fj,B, as obtained by fitting to the result for TJ^Ms) Htl and 
demanding that Q = 0AB and S = 0. 



The overall grand-canonical reference environment is 
determined as follows. Using the baryon chemical po- 
tential hb as a control parameter, we obtain the freeze- 
out t emp erature Th from the fit to the data obtained in 
Ref. [3, yielding T h (» B ) » T [l - (fi B /m N ) 5 / 2 } with 
T = 170 MeV. Subsequently, we perform a grand- 
canonical iteration to determine those values of /xq and 
fis that ensure Q = aB and S = 0, where a = 0.4 which 
is representative of Z/A for gold. The resulting values are 
shown in Fig. ^ as functions of [ib- As \ib is increased, 
baryons become favored over antibaryons, so there will a 
bias of hyperons over antihyperons. To counterbalance 
the associated net negative strangeness (recall that Sa is 
minus one), /ig must take on a positive value to ensure 
a compensating excess of kaons over antikaons, leading 
to an approximate proportionality between ^5 and fiB- 
This balancing of the strangeness in turn leads to a small 
bias towards positive net charge and fj.Q must therefore 
be correspondingly negative. While the resulting value 
is also approximately proportional to fis, the absolute 
value is rather small and one might for most purposes 
simply take (iq to be zero. 

In our spinodal scenario the value of So is sampled at 
the plasma stage and then canonically conserved though 
the hadronization process. Each plasma blob has the vol- 
ume V q (usually taken to be 50 fm 3 ) and its strangeness 
content is determined by sampling v and v at the plasma 
temperature T q . which is either taken to be equal to Tq 
(= 170 MeV) or to equal the hadronic freezeout temper- 
ature T/j(/is), which should provide approximate upper 
and lower bounds. At fiB = 0, where T q = 170 MeV, 
the width of the strangeness distribution in such a blob 
is as = 5.74. If we take T q = Th, the distribution will 
grow narrower as [Ib is increased and at [Ib = 400 MeV, 
where the temperature has decreased to 150 MeV, we 
have as = 4.65. 

The blob then expands until the freezeout volume 
Vh = "xYq is reached. (Guided by rough approximations 



FIG. 2: The average K multiplicities as functions of the 
baryon chemical potential /is in three scenarios: 1) the stan- 
dard grand-canonical treatment where the hadronic freeze-out 
occurs in the volume Vh at the temperature T/^/xb), 2) our 
spinodal treatment in which the blob strangeness So is sam- 
pled in the plasma volume V q at the temperature T q = To 
and then conserved through the hadronic freeze-out, and 3) 
the restricted canonical treatment admitting only So = 0. 



to the equation of state, we usually employ an expansion 
factor of x — 3, i.e. Vh = 150 fm 3 .) At this point, the 
blob has transformed itself into an assembly of hadrons 
whose abundances are assumed to be governed by the 
canonical distribution (|15fl associated with the particular 
value of So, and the modified chemical potentials H B (So) 
and /Ltg(So) (see Eq. I|iup) that have been adjusted for 
each particular blob strangeness So to ensure matching 
of the corresponding biased baryon and charge contents, 
as explained earlier. 



A. Kaon multiplicity distributions 

We first discuss the results for the charged kaon mul- 
tiplicity distributions. Figure |21 shows the resulting aver- 
age multiplicities for three different scenarios. The 
first scenario is the usual global grand-canonical treat- 
ment (see above), while the second is our spinodal sce- 
nario in which the blob strangeness is determined at the 
plasma stage and then kept fixed during the hadroniza- 
tion. Both of these scenarios consider all possible values 
of S*o but while the distribution of this quantity is de- 
termined at the hadronic freeze-out in the former, it is 
determined already in the plasma in the latter. In the 
third "restricted canonical" scenario the blob strangeness 
is always required to vanish, So = 0. 

The various scenarios lead to rather similar results. 
The K~ multiplicity decreases steadily as a result of 
the strangeness balancing explained above while the K + 
yield initially increases, for the same reason, but then, 
as the freeze-out temperature Th{^s) begins to decrease 
noticeably, the overall hadron production decreases, thus 
yielding a decreasing behavior of the K + curve. Gener- 
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FIG. 3: The variance in the K + multiplicity as functions of 
the baryon chemical potential /i_g for the various scenarios 
considered in Fig. using either To or Th for the plasma 
temperature T q in the spinodal scenario. 

ally, the average multiplicities are affected very little by 
the blob formation, relative to the grand-canonical treat- 
ment, and the spinodal results are therefore omitted for 
K~ while only the results for T q = Tq are shown for 
K + . These exhibit a progressively increasing enhance- 
ment that amount to 2% at {is = 400 MeV. By con- 
trast, for the restricted scenario, where only blobs having 
Sq = are admitted, there is a reduction in the averages 
by a few per cent for all hb- 

The picture changes when the multiplicity fluctuations 
arc considered, as illustrated in Fig. |3J where the corre- 
sponding multiplicity variances are shown for the posi- 
tive kaons. While the overall behavior is qualitatively 
similar to the behavior of the averages, there are sev- 
eral important differences. First, in the restricted sce- 
nario (where only Sq = is included) the suppression 
of the variance is significantly larger than was the case 
for the average, amounting here to 8-10%. Furthermore, 
at the larger values of (is, where the net baryon density 
becomes significant, the grand-canonical variance suffers 
more from the decreasing temperature than the spinodal 
variance. This important divergence is a result of the fact 
that the larger average baryon number implies a corre- 
spondingly larger baryon-number variance as well and 
therefore also a larger variance in the strangeness (since 
strangeness in the plasma is carried exclusively by quarks 
and antiquarks which also carry baryon number). Most 
importantly, there is a significant dependence on the em- 
ployed plasma temperature T q . In order to illustrate this 
effect, we show results for two extreme values, namely 
T q = T = 170 MeV and T q = T h (fj, B ), and for these the 
difference amounts to ten per cent at hb = 400 MeV. 

B. Multiplicity ratios 

From the experimental perspective, it is more conve- 
nient to consider multiplicity ratios, and we therefore 
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FIG. 4: The average value of the K/n ratio for either posi- 
tive (increasing) or negative charges (decreasing) for the three 
scenarios considered in Figs. 121. il 
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FIG. 5: The normalized variance in the K/n ratio for either 
positive (increasing) or negative charges (decreasing) for the 
three scenarios considered in Figs. 12131 The variances have 
been multiplied by the corresponding average pion multiplic- 
ity in order to make the result scale invariant. 

show in Fig. 0] the average K + /n + and K~/n~ ratios, 
for the same three scenarios. When fi B is positive there 
is a preference for K + over K~ and hence K + /n + will 
increase while K~ /n~ decreases. This behavior is prac- 
tically linear since the suppression from the decreasing 
temperature affects all hadrons species. Since there is a 
(small) tendency for the n and K multiplicities to vary 
in concert, the difference between the various scenarios 
is reduced. In particular, there is hardly any difference 
to be seen for the K~ /n~ ratio. 

We now turn to the corresponding variances. Since 
the variance of the K/n ratio decreases in inverse pro- 
portion to the size of the system, it is convenient to 
multiply by the mean pion multiplicity and thus obtain 
a result that approaches a constant for large volumes, 
(7r ± )(T 2 (/\ ± /7r ± ). The resulting results for the positively 
charged hadrons are shown in Fig. [SJ They are quali- 
tatively similar to those for (K /tt ) (Fig. 0J. But al- 
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plasma volume V q . Since both averages and variances 
are proportional to volume, a smaller V q shrinks the dis- 
tribution of the blob strangeness So (so it has smaller 
mean and variance). Consequently, the resulting values 
(K + /ir + ) and (tt + )(j 2 (K + / i tt ± ) are decreasing functions 
of x- However, this dependence is not dramatic: a dou- 
bling of x from 2 to 4 reduces the average ratio by less 
than 2% and the normalized variance by less than 10%. 

It is important to recognize that there are two oppos- 
ing effects: One is the basic fact there are more degrees 
of freedom in the deconfined plasma phase than in the 
confined hadron gas, which enhances the fluctuations. 
(At T = T and fi B = we have ( s « 0.33 V q while 
£s=±i « 0.081 Vh, so the effective degeneracy of the 
plasma is approximately four times larger than that of 
the hadron gas, which would then be compensated with 
an expansion factor of x ~ 4, as indeed borne out by 
the calculation.) However, as x is increased from unity, 
this advantage is being steadily offset by the ever larger 
volume of the freeze-out configuration, Vh = xVq- The 
crossover happens to occur at x ~ 3.3, a value rather 
near our adopted estimate, x = 3. Should it turn out 
that there is less change in volume from the formation 
of the plasma blob to the hadronic freezeout, then the 
relative effect of the strangeness trapping will be consid- 
erably larger. For example, if x were only ten per cent 
smaller, the effect would be about twice as large. 



FIG. 6: The average of the K + /tt + ratio (top) and its variance 
multiplied by the mean 7r + multiplicity (bottom) as obtained 
at fiB — 300 MeV for the scenarios considered in Fig. El shown 
as a function of the volume expansion factor \ = Vh/V q . 



though variances in the ratios are less sensitive to the 
specific scenario than the kaon variances themselves (Fig. 
|3J), the differences are still clearly brought out. Further- 
more, as expected from Fig. [31 an increase of the plasma 
temperature T q increases the resulting values. The fluc- 
tuations in the K + /ir + ratio may thus offer a suitable 
observable that is sensitive to a dumping-induced trap- 
ping of strangeness in the expanding matter prior to the 
hadronization. 



C. Dependence on the expansion factor 

The above results have been obtained for a given ex- 
pansion factor x — Vh/Vq = 3. This value should be 
taken only as a rough approximation to what might ac- 
tually happen and since the results are sensitive to this 
parameter it is of interest to consider also other de- 
grees of expansion. This aspect is illustrated in Fig. 
where the average of the K + /ir + ratio and its variance 
(normalized by (tt + )) are shown as a function of x f° r 
Hb = 300 MeV. Since we keep the hadronic freezeout 
volume equal to Vh = 150 fm 3 to facilitate the com- 
parisons, a larger expansion factor x implies a smaller 



D. Effect of source mixing 

The above studies have considered only the hadrons 
resulting from a single blob which, presumably, popu- 
late a certain limited kincmatical region centered around 
the velocity of the original blob |12|. However, even if 
a spinodal decomposition into kinematically separated 
blobs were to occur, it would not be experimentally fea- 
sible to restrict the measurement to include only those 
hadrons resulting from a single blob. Rather, one should 
generally expect that a given detection acceptance will 
admit hadrons originating from more than one single 
blob. Thus one needs to address the fact that a mix- 
ing of hadrons from different sources will degrade the 
strangeness-trapping signal. 

In order to elucidate this practically important feature, 
we calculate the production from several different blobs 
and combine the resulting hadrons into a single "event" 
before performing the analysis. The resulting variance in 
the K/tt ratio (multiplied by the mean 7r + multiplicity) 
is shown in Fig. [7| as a function of the number of blobs 
whose products have been combined. While the value 
drops by about a factor of two when going from a single 
blobs to two combined blobs, it subsequently stabilizes 
and quickly approaches a constant as ever more blobs are 
combined. Furthermore, the relative increase when go- 
ing from the standard grand-canonical scenario one of our 
spinodal scenarios remains nearly unaffected by the num- 
bers of blobs combined in the analysis. This result indi- 
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FIG. 7: The effect of combining the hadrons from several 
separate plasma blobs before extracting the variance of the 
K + /% + ratio, for the three different treatments. (For these 
results we used an expansion ratio of \ = 2.) 
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FIG. 8: The dependence of the coefficients n avc , <r var , a cov 
on the total number of blobs that are subject to the global 
strangeness constraint, Stot =0, for K + in a scenario having 
fj,B = 300 MeV, Tq = T , and x = 2, as in Fig. 



cates that the suggested signal of the strangeness trap- 
ping is robust against the inevitable source mixing and, 
consequently, it may in fact be practically observable. 



E. Global strangeness conservation 

We finally analyze in some detail the role played by 
the overall conservation of strangeness in each event. 

In the above studies, we have treated the strangeness in 
a given plasma blob as a grand-canonical variable, which 
is well justified when the blob volume forms only a small 
part of the total system. In order to investigate the qual- 
ity of this approximation, we consider the effect of con- 
straining the combined strangeness of ./V individual blobs 
to zero. The corresponding canonical partition function 
for the combined system is then 



r (i-JV) 
-tj 



N 

n 



Si 



Ss!+...+s N ,o , (26) 



(i) 

where Z s is the canonical partition function for blob 
i having the strangeness Si. For simplicity, we shall 
assume that all TV blobs are similar (as we have done 
above). First, to obtain the correlated distribution of the 
blob strangenesses (Si, . . . , Sn), we consider the problem 
at the plasma level where Zq — Iq(2N( s ). We then find 



iV 



Po(S±, . . . , Sn) 



J5i + ...+Sjy,0 

/o(2TVC s 



] [is, (2Cs)] • (27) 



as described in more detail in Appendix IO 

In order to bring out the effect of the global strangeness 
constraint, we consider the multiplicity of a particular 
strange hadron k resulting from the combined observa- 
tion of some of the blobs, say those labeled 1, . . . , TV', 



i.e. n* 



where tin is the multiplicity 



contributed by the blob i. If all the blobs are similar, the 
average multiplicity is of the form 



<n K ) = TV'nr (N) 



(28) 



where n^ ve (TV) is the average multiplicity arising from 
any single blob. Furthermore, the variance in the total 
multiplicity n K has the following form, 



N'[ar(N) 



TV' - 1 
TV- 1 



(29) 



where cr™ r (TV) is the variance in multiplicity from any 
single source and cr£ ov (TV) is the covariance between the 
multiplicity from any one source and the combined mul- 
tiplicity from the all the other TV — 1 sources. The overall 
restriction of the total strangeness to zero reduces the 
individual multiplicities {n! }, leading to smaller val- 
ues of both n^ ve and u™ r . Furthermore, a higher-than- 
average strangeness in one blob introduces a bias towards 
lower-than-average values in the others, in that particu- 
lar "event" , thus producing an anticorrelation among the 
individual strangeness values {S n }. This in turn leads to 
negative covariances, a c ° v < 0. 

To illustrate these effects, we consider the production 
of positive kaons and show in Fig.[5]the dependence of the 
coefficients n avc , er var , and cr cov on TV, the total number of 
blobs that are subject to the global constraint. (We may 
here employ either analytical recursion relations or statis- 
tical simulation, as discussed in Appendix IU1 ) For TV=1 
the strangeness of each blob must vanish and we obtain 
the restricted canonical scenario considered earlier. In 
the opposite extreme, TV — > oo, the global constraint be- 
comes ineffective and the standard grand-canonical sce- 
nario is recovered. To a good approximation, the devia- 
tions of the coefficients from their grand-canonical values 
are inversely proportional to TV, for example, 



r (TV) 



r (oo) 



1 

TV 



:^ ar (i) 



ar(oo)]. (30) 
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FIG. 9: The ratio of the variance to the mean for the K + 
multiplicity resulting from a combined system of TV blobs of 
which TV' are being observed, plotted as a function of TV' /TV. 



These results were obtained for jj,B = 
with either x = 2 (upper curve) or x 



FIG. 10: The variance of the K + /ir + ratio for a single blob, 
multiplied by the mean n + multiplicity, when the combined 
net strangeness of TV blobs is required to vanish. These results 



300 MeV and T q = T wer e obtained for /i B = 300 MeV and T q = To with either 



3 (lower curve). 



X — 2 (upper curve) or x = 3 (lower curve), as in Fig. El 



One may then readily judge the effect for a given TV on the 
basis of the change in going from TV = oo to TV= 1. For 
the average number of K + emitted by each blob, n avo , 
this drop is only about 4.5% (2.5%), while it amounts 
to 25% (16%) for the corresponding variance, er var (the 
numbers in paranthesis refer to x = 3) ■ This behavior may 
be contrasted with the approximate constancy of cr var + 
cr cov , the variance of the total K + multiplicity divided 
by TV, which decreases by only 4% (2%) when going from 
TV ^ 1 to TV = 1 . These results illustrate the fact that the 
variance in the total multiplicity is always smaller than 
the corresponding mean when the system is subject to a 
canonical constraint. 

Figure [U] illustrates the dependence of the canonical 
effect on the ratio TV'/TV by displaying the ratio of the 
variance to the mean for the K + multiplicity resulting 
from a combined system of TV blobs of which TV' are being 
observed. For any given value of the expansion ratio \ = 
Vh/Vg, which governs the effective statistical weight of 
the plasma relative to that of the hadron gas, the change 
from the limit where only a small fraction of the combined 
system is being observed to the situation when all the 
kaons are collected is seen to follow a universal curve 
that is linear in TV'/TV, to a very good aproximation. 

The corresponding effect on the fluctuations of the 
K/tt ratio is illustrated in Fig.^] It shows the value of 
(tt + )(j 2 (K + /tt + ) obtained for the hadrons emitted from 
a single blob when TV such blobs satisfy a combined 
canoncial constraint on their total strangeness contents, 
Stot = 0. The canonical constraint is obviously most 
effective in the extreme case of a single blob, TV = 1, 
where it is identical to the restricted scenario in which 
only So = is allowed for each individual blob (bot- 
tom curve in Fig. EJ) . As more blobs share the burden of 
the constraint, its effect rapidly diminishes and the re- 
sults approach those of our standard spinodal scenario, 
where the strangeness of each individual plasma blob is 



determined grand canonically (top curve in Fig.0); this 
approach occurs approximately as ~ 1/TV. Thus already 
at TV w 5 the effect is hardly significant. Consequently, if 
a given blob represents less than 20% of the total system, 
the vanishing of the overall strangeness should have no 
noticeable effect on the observables considered here. 



V. CONCLUDING DISCUSSION 

Most observables considered so far in high-energy colli- 
sion experiments behave rather smoothly as functions of 
the control parameters, such as bombarding energy and 
centrality. However, preliminary data analysis of fixed- 
target experiments at the SPS by the NA49 collaboration 
have revealed two intriguing exceptions. First, there 
appears to be a significant enhancement of the K/tt ratio 
at beam energies of 20 — 30 ^4GcV. Second, in the same 
energy range, the fluctuation of this ratio (expressed rela- 
tive to mixed events) is strongly enhanced |2ll | . Both en- 
hancements are localized within about 10 AGeV in beam 
energy. Though otherwise rather successful in describ- 
ing the hadron yields, statistical models cannot repro- 
duce the observed energy dependence of these enhance- 
ments • Nor does a hadron gas and its resonance- 
induced correlations account for the strong enhancement 
in the ratio fluctuations [2^, 01 ■ O ne might then won- 
der whether these strong fluctuations might result from 
the enhanced fluctuations associated with a second-order 
phase transition However, if this were the case, one 
would expect particularly strong fluctuations in the pi- 
ons. However, while the fluctuations of the K/n ratio are 
observed to be enhanced, those of the p/ir ratio follow the 
expectations from transport models |20| . 

If the statistical models 0] indeed provide reason- 
able estimates for the relation between beam energy 
and the thermodynamic parameters characterizing the 
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chemical freeze-out, i.e. the temperature and the chemi- 
cal potentials, then the anomalous behavior reported by 
NA49 would be consistent with a first order-phase tran- 
sition having occurred prior to the chemical freeze-out, 
at a somewhat higher temperature and about the same 
chemical potential. In order to evaluate the plausibil- 
ity of such a speculation, we have in this paper stud- 
ied the effects on yields and fluctuations if the bulk of 
the system indeed undergoes a spinodal-like break-up as 
it hadronizes. The key physical effect is then that the 
amount of strangeness residing within the plasma that 
forms a given blob remains effectively trapped, thus im- 
posing a canoncial strangeness constraint on its conver- 
sion into hadrons. 

Considering idealized scenarios in which an equili- 
brated plasma forms a number of separate blobs that sub- 
sequently expand and hadronizc independently, we have 
developed the relevant statistical tools. These results 
may be of general applicability. With a view towards 
the specific observations at the SPS, we have studied the 
production of pions and kaons in such a scenario rela- 
tive to the standard picture in which statistical equilib- 
rium is maintained through freeze-out. While the average 
hadron yields are essentially unaffected by the breakup, 
the fluctuations in the strange-particle multiplicities are 
significantly enhanced, especially at larger values of the 
chemical potential. From the experimental perspective, 
the yield ratios are particularly interesting and we have 
especially studied the K/ir ratio. Depending on the de- 
gree of expansion during the stage of strangeness trap- 
ping, the fluctuations in the K + /tt + ratio can be en- 
hanced significantly (of the order of 10%) relative to 
the standard scenario where global equilibrium is main- 
tained. 

Thus, the present (rather idealized) studies suggest 
that a spinodal decomposition might indeed lead to en- 
hancements of the magnitude observed by NA49. How- 
ever, before such a conclusion could be made with con- 
fidence, further studies would be required. In particu- 
lar, both strong resonance decays and weak decays must 
be taken into account. Moreover, the enhancement of 
the K/ir fluctuations should be correlated with other 
expected effects, such as A^-body kinematic correlations 
[12] ■ Of particular interest are more precise calculations 
of the equation of state in the relevant baryon-rich envi- 
ronments, to better ascertain the location of the expected 
phase coexistence region, as well as refined calculations 
of the collision dynamics to determine whether the spin- 
odal region is in fact likely to be encountered and, if so, 
to what degree the phase decomposition may actually 
develop. We hope that the present findings will provide 
added incentive for these challenging undertakings. 



APPENDIX A: CANONICAL FERMIONS 

We discuss here the canonical treatment of fermions, 
in which the total number of quanta present is kept fixed 



to a given value Nq (the number of "particles). The ap- 
propriate partition function is then 



ZnM = ] 



E 



-fluid 



E 

{™i}« 



-/3E{m} 



SiY^rii-No) (Al) 

i 

, E{m} = ^2 riiSi . 



Here Ui is the number of quanta present in the par- 
ticular "single-particle" state i of energy Further- 
more, {ni}No denotes the subset of all configurations {rii} 
whose total particle number is constrained to be equal 
to the specified value No, N{rii} = J2i n i = -^o- The 
fermionic restriction on n to be at most one complicates 
the evaluation of the partition function. However, it is 
possible to derive the following recursion relation, 

1 N 

ZnW = -^R-^j^zfM) , (A2) 



with Zq(/3) = 1. N can be readily sampled numerically 
on the basis of P(N) = Z^/Z, where Z F = J2n>o^n 
is the grand-canonical partition function. 

One way to derive the above recursion relation reorga- 
nizes the basic expression for Zn, 



J N 



E 



U<«<"'<W 



- y 

N\ ^ 



3 a il +a i3 +-+a ilt ( A g) 
e a si +a l2 +---+a ijv ^ 



il#i2^---#ijv 

where = — (3ei and no indices may be equal in the 
second sum. With £„ = Zi(n/3), the first few terms are 



\ [Z.iP) 2 - Z 1 (2/3)] = \ [ZxCi - Z0C2] , (A5) 



3! 



3! 



i[£2Cl-£lC2+£oC3] , 



e ai+aj+a k 

ijk 



(A6) 



The first term is the classical limit, Zn = Z± /N\ H — ■ 
Another approach starts from the expansion of InZ, 



InZ 



= EMi 

i 

= E 



EE(- 

i n>0 



n 



n-l 



■(n — Cl — 5C2 + qC3 



n>0 



then exponentiates InZ = [Ci — 5C2 + 5C3 + 



y-i 



N>0 



l + ]((l+i[. ]((! 



,(A7) 

(A8) 
)))). 



11 



and finally reorganizes the series in powers of e^ £ , 

Z = Y, z n = 1 + Ci + 5K1-C2] (A9) 
-° z 2 

+i[i(C 1 2 -C 2 Ri-CiC2 + C3] + ... • 

a. Degeneracy and volume. The above derivations 
have been made for non-generate systems, g = 1. In the 
general case, when g > 1, the situation is more compli- 
cated since a given energy level may accommodate up 
to g quanta. There is then no simple relation between 
Z[g = 1] and Z' = Z[g > 1]. However, using the above 
expansion (|A7|) of In Z in terms of £„ we find 

InZ' = . 9 ^ln[l + e-^] = ... = , 

(A10) 

where (? n = gQ n = gZ\ (n(3) . The partition function may 
thus be obtained by replacing ( n by (' n = gC, n in the above 
procedure, 

z' = Y, z 'n , z'n = Zn{Cu - C = gCn] . (Ail) 

Consequently, the recursion relation takes the form, 
1 N 

^ = ¥EH n X-nC,2o = l, (A12) 

71=1 

with = = gJ2t e~™' 9c * , and -2^- becomes an N th 
order polynomial in the degeneracy g, 

Z' = 1 , Z[ = gd , Z^ = \g 2 d - \g(, 2 ,(A13) 

z 3 = y 3 ci - WC1C2 + y< 3 , ■ ■ ■ ■ (Au) 

We note that the dependence on the volume V is similar 
to the dependence on the degeneracy <?, since both enter 
as overall factors in the elementary partition functions, 

Cn ~ gV. 

b. Multiplicity distribution. The mean multiplicity 
is enhanced/reduced for bosons/fermions, 

Z x 2tt 2 h 3 c 3 ^ ' n 3 y T 

= ^(±)"- 1 C„ = Ci±C2 + --- , (A15) 

n>0 

and so is the corresponding multiplicity variance, 

2 g VT 3 ^ j 1 - nm 

a N = TTTFT — K 2{^r) (Alb) 

27T h 6 c 3 *—t n z T 

= z ^(±)"- 1 r<„=Ci±2C 2 +3C 3 ±--- , 

n>0 

and both are strictly proportional to both the volume 
V and the degeneracy g, since £ n ~ gV , as just noted 
above. 



c. Chemical potentials. The above results apply in 
the absence of a chemical potential. The presence of a 
chemical potential effectively replaces the energy ej by a 
shifted value — fi and consequently all the manipula- 
tions can be carried through as above. So, with 

( n (a,f3) = UP)o- an , U(3) = 2iM =J2e~ n ^ , 

(A17) 

we see that all the terms in Zn contain the same power 
of the fugacity e~ a and we find 

Z N (a,0) = Z N {C°(n(3)}c- aN . (A18) 

However, the multiplicity moments have a complicated a 
dependence, 

(AO = Ci e" Q ± C2 c- 2a + d e~ 3a ± . . . , (A19) 

al = C°ie""±2C2e- 2Q +3C3e- 3a ±-- - . (A20) 

Furthermore, the corresponding expressions for the asso- 
ciated antiparticle are 

(AO = Ci e +a ± C2 c +2a + d e +3Q ± . . . , (A21) 
0-% = (ic +a ±2 ( 2 c +2a + 3 ( 3 e +3a ± • • • . (A22) 

In the case of a plasma blob, which has no strangeness 
bias, the distributions of N and N must be identical. 
Consequently, in that situation, we must have a = 0. 
Thus, if the s and s quarks are embedded in an en- 
vironment where fis (and/or /iq) have finite values, a 
strangeness symmetric distribution can be established by 
adjusting /is appropriately: /is = — A*q) (since 

B s = §, Qs = and S s = -§). 

It may seem odd to introduce a chemical potential 
within the context of a canonical treatment, but a given 
scenario may well require a canonical treatment with re- 
spect to one attribute (e.g. strangeness), while a grand- 
canonical treatment suffices with respect to another (e.g. 
baryon number). Of course, as the above analysis brings 
out, the consideration of several attributes is relevant 
only when the system contains several species that com- 
bine the attributes differently. (The inherent correlation 
between baryon number and strangeness in the quark- 
gluon plasma was recently proposed as a useful diagnostic 
for strongly interacting matter [25|.) 

APPENDIX B: GENERAL STATISTICAL 
TREATMENT OF STRANGE HADRONS 

We derive here the expressions needed for the general 
(classical) statistical treatment of a gas of hadrons char- 
acterized by a temperature T and chemical potentials for 
baryons and electric charge, [Lb and [1q, with its total 
strangeness So being kept fixed. The strategy will be to 
group the strange hadrons species {k} together according 
to their strangeness S K and then build up the complete 
partition function by pairwise inclusion of classes having 
opposite value of their strangeness. 
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1. Generic multiplicities 

The one-particle partition function for a given strange 
hadronic specie k is given by 

MT,n B ,n Q ) = £ K (T) ^bB^QJ/t , (Bl) 

where its value for vanishing chemical potentials is 



g K VT 3 fm K c 2 ^ 2 - 2 



2tt 2 h 3 c 3 V T 



A" 2 



(B2) 



The effective one-particle partition function for a class of 
hadrons having a common strangeness S is then 

Cs(T,fiB,fiQ) = 22^s K ,sC K (T,fJ,B,fJ.Q) (B3) 

K 

= X> re ,s C,(T)e^ B -+"««-)/ T , 



which generally differs from C-S- The number of hadrons 
having the strangeness 5 in a given system is denoted by 
Ns and the probability that such a hadron belongs to 
the particular specie k is given by Ct/Cs- The associated 
first and second multiplicity moments are then (n K >) = 
(Ns'Krt/Cs' and (n K .n K ») = (N S ,N S »)( K >( K » /Cs<s" ■ 
We may therefore concentrate on finding the generic mul- 
tiplicities {N$}. 

We first combine two conjugate classes {+S 1 } and {— S} 
to form the class {±>S}. The partition function for the 
resulting combined system of hadrons having S K = ±S 
is then 



A±s} 

'So 



E 



+S 



! TVT^I °( w +s-^-s)S,So 



k So 



W2Co) , 



(B4) 



where ./V±5 > denotes the number of hadrons having 
the strangeness S K — ±S and £q = C+sC-s- We note that 

{±5} 

Zg Q and the corresponding multiplicity moments van- 
ish unless So is a multiple of S, i.e. So — 0, ±5, ±2S, .... 
The factorial multiplicity moments can also readily be 
obtained, 



(N +S (N +S 

/-m ( C—S 

VC+s 



-1) 



. . (N+s - m 

1 So -mS ( 2 Co) 
^5 (2Co) 



D)£ s} 



Furthermore, the mixed multiplicity moments can be ob- 
tained by use of the constraint (N + s — N-s)S = So- 

These relations provide a complete treatment of one 
pair of conjugate classes. Imagine now that we have 
thus obtained the partition functions Zg and the cor- 
responding multiplicity moments (N± s ) s for \S\ = 



1,2,3. The classes {±5} may then be combined recur- 
sively. Thus, combining first {±1} with {±2}, we find the 
partition function for the combined ensemble {±1, ±2}, 



r{±l,±2} 



V- „{±1} -{±2} 



""So — So, S2 



(B6) 



S 2 =0,±2, 



where So is the specified total strangeness. It follows that 
if we have a self-conjugate system (i.e. for each hadron 
specie k included, the corresponding antispecic R is also 
included) with strangeness S K = ±1,±2 and whose total 
strangeness is So, then the probability that those with 
S K = ±1 have a combined strangeness of S' is given by 

p{±l,±2} {s{±iy= gl) = Z {±1} Z {±^ /Z {±1, ±2}; m 

which vanishes unless So — S' is even. Consequently, 
after the classes {±1} and {±2} have been combined, 
the multiplicity moments for hadrons with S K = ±1 are 



(NF=±i) { st h±2} 



E r 

S' — Sq ,Sq i2 , 



(4t h±2} Y (B8) 

^{±1} ^{±2} /Afm x{±l} 
■°S' ^So-S' \ JV S=±1/S' ' 



while those for hadrons with S K = ±2 are 



v{±l,±2} 



±2 /So 

E 

S 2 =0,±2,±4, 



= z 



,{±1,±2} 
'So 

,{±1} ^{±2} 



(B9) 



x{±2} 
S=±2/S 2 



Similar expressions hold for the mixed moments, e.g. 

<Au;v +1 ><f ±2 > = (4 ±1 ' ±2} )" ( B1 °) 

x £ z * 1} z t-i <^+i>r } • 

S' — So , So i2 , . . . 

It is straightforward to verify the following sum rule 
expressing strangeness conservation, 



E w 

S=±l.±2 



,{±1,±2} Q _ 



S/S 



S — So 



(Bll) 



We note that the above recursion scheme holds even 
if there are no hadrons with S = ±2. In that case, 
the corresponding effective partition functions vanish, 
(,±2 = 0, and, as noted earlier, the combined partition 
function is unity for Sq = and vanishes otherwise, 



'{±2} 
'So 



5s 0i o- As a consequence, the combined par- 
(B5) tition function remains unchanged by the incorporation 
of S = ±2, zj^ 1 '^ 2 ^ = Zq^ 1 ^, and the multiplicity mo- 



So So 

mcnts remain unchanged as well. 

Further conjugate strangeness classes may be incorpo- 
rated analogously. Thus, the inclusion of S = ±3 yields 
the following partition function for {±1,±2,±3}, 



; {±1,±2,±3} = Z {±1,±2} 2 {±3} 

"''S'o / J Sq> — S3 S3 

S 3 =Q,±3, 



(B12) 
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and the generic multiplicity moments arc given by 



/Arm \{±1,±2,±3} 
\^S=±l)s - [f'So 

x £ z£% (N F=±1 )^ , 

S'=S ,So±3,... 

x £ ^ 1,±2} 4^ (^ ±2 )l ±1 ' ±2} , 



r {±l,±2,±3} 



(B13) 



S'=So,S ±3,... 



(B15) 



>< e 4 3 ±3} (^= ±3 >r } 

S 3 =0,±3,±6,... 



and the sum rule becomes 

e ™r ,±2i±3} 

S=±l,±2,±3 



5 = So 



(B16) 



Recursion expressions for the correlations between the 
generic multiplicities can be obtained in a similar man- 
ner. For example, the correlations between iV_i and iV +1 
follow from the corresponding mixed moment, 



Y: 4^' ±2} 4S'^-^4^ {±1>±2} 



S ; — So , So i3, 



This procedure readily generalizes to the combination of 
any number of self-conjugate classes. 



The second multiplicity moments are more compli- 
cated to obtain, because they receive contributions both 
from the correlated internal multiplicity fluctuations 
within the separate S classes characterized by each par- 
ticular generic multiplicity set {-/Vg} and from the corre- 
lated fluctuations of these generic multiplicities. It is con- 
venient to introduce the correlation coefficient o K ' K " = 
$S'S" Pk' {$k'k" ~ Pk") which expresses the degree of cor- 
relation between two particles in the same strangeness 
class. Thus it vanishes if S" ^ S" . When S' = S" 
we have a K ' K " = —p K >p K " when the two species differ, 
k' =/= k" , while <j k i k i> = p K >q K ' when k! = k" , where 
p K i = ( K i I (s' is the probability that a hadron of the class 
S' belongs to the specie k' and q K > = 1 — p K < is its com- 
plement. We then find 

(n K >n K// ) So = P K 'P K "(N S 'N S "}s + ^S'S"(Ns'}s ^ k ,' k " ■ 

(B20) 

The multiplicity covariances then involve both the intra- 
class correlations un'N" arid the inter-class correlations 

CTri'n" = Ss'S"{N')s ^K'K"i 



,S 



(n K >n K „) Sa - (n K >} So (n K »}s (B21) 
Pk>Pk"Vn>n" + Ss'S" (N') So cr K ' K " . 



The above expressions allow us to evaluate the average 
multiplicities of individual hadron species as well as the 
associated (co)variances, for any given value of the fixed 
strangeness of the blob, Sq. 



2. Individual hadron species 

The above treatment allows us to determine the canon- 
ical moments of the generic hadron multiplicities {iVs} 
in any blob with a specified strangeness So. We now con- 
sider the multiplicities of the individual hadron species. 

For a given value of Sq, the mean number of a partic- 
ular species k' is given by 



(n K >) 



So 



1 



n 



E 

n K >0 



-S ) , 
(B18) 

and the correlation between the multiplicities of any two 
species k' and k" is given by a similar expression with 
n K ' replaced by n K 'n K ». In order to use these expres- 
sions with the above expressions for the generic multi- 
plicities {iVc?}, which no longer contain the individual 
species multiplicities {n K }, we note that the probability 
that a generic hadron of strangeness S — S K is of a par- 
ticular specie k is given by p K = Ck/Cs k - Consequently, 



/So 



p K >(N. 



S K ,)S 



(B19) 



APPENDIX C: GLOBAL CANONICAL 
TREATMENT 



We consider here the more complicated situation where 
N individual blobs arc subject to a global canonical con- 
straint on their combined strangeness, Stot = ^2 n ^n- 
We specialize to the relevant case of Stot = 0, but the 
treatment can readily be adapted to any value. 

Assuming, as we have throughout, that all the blobs 
are similar, the joint probability for finding the combined 
system with v n s quarks and v n s antiquarks in the blob 
n is then given by 



P(v 1 ,v 1 \...\v N v N ) = 
I (2NCs) vi\ vi 



(CI) 



J V\A \-VN,Vl"- + VN ■ 



This expression can be used as a basis for a di- 
rect Metropolis sampling of the individual multiplicities 
\y n , starting for example from {0, 0}. However, this 
may not be optimally efficient, since in fact we only need 
to know the distribution of the differences S n = v n — v n . 
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which is given by 



P(S!,...,S N ) 

N 

n 



E S Vn-<>n 



P(i/ 1 ,i/ 1 ;...;v n l> n ) 



j Si + ---+Sn,0 



N 

n 

(JSi + .-.+Sjv.O 



E 



I [J S „(2C.)] 



(C2) 



The corresponding Metropolis sampling could then start 
from {S n } = {0}, for example, and repeatedly exchange 
one unit of strangeness between two selected subsystems, 



Si 



S : 



Si =b 1 and Sn 



Sj =F 1, with the 



corresponding ratio of weights being given in terms of 
ratios of Bessel functions of neighboring order, W /W = 

(Is t ±i(2Cs)/IsMsWs jT i(2C s )/IsA^s))- 

A much simpler procedure is based on the fact that 
the entire set of configurations {v n , v n } can be organized 
according to the total number of s quarks present, M 
(which equals the total number of s when Stot = 0). The 
class having M = contains only the empty configura- 
tion {0,0}, while the class having M = 1 has N 2 mem- 
bers, and so on. Thus, class M has (N M /Ml) 2 members, 
each of which has the relative weight ( 2 M , and it is read- 
ily checked that the sum of weights is Io(2N( s ), the total 
partition fucntion. The expected number of s quarks is 
N( s and the most likely value of M is [iV£ s ]. It is easy 
to sample M by a Metropolis procedure based on the 
weight ratios for adjacent values of M, Wm+i/Wm = 
(N( s /M + l) 2 and W M -i/W M = (M/N( s ) 2 . 

Alternatively, M could be sampled directly from its an- 
alytical distribution P(M) = ( 2M ' (N M /Ml) 2 / I (2N( S ) 
in a standard manner. Once M has been selected, it is 
straightforward to distribute the M quarks and M an- 
tiquarks randomly among the N subsystems and thus 
obtain their strangenesses as S„ = u„ — u„. 

The above discussion brings out the fact that the total 
plasma partition function can be calculated by either dis- 
tributing v n quarks and v n antiquarks in each subsystem 
n or distributing M quarks and M antiquarks anywhere 
within the combined system, 



N 

n 



E 



v n \ v n l 



fai+. 



(C3) 



2^ — ~' 77T- °M,M = loKZNQs) 



MM 



Ml Ml 



In any case, once the individual blob strangenesses 
{S n } have been selected, the various hadron multiplic- 
ities can be sampled as in the standard case discussed in 
Appendix It is thus possible, in this manner, to make 
a complete statistical sampling of the combined system 
and then perform any analysis of interest. 



However, such a full similation is not always neces- 
sary. In particular, as is often the case, when the quan- 
tities of interest can be expressed in terms of first and 
second moments of the multiplicity distributions of spec- 
ified hadron species, it is possible to employ recursion re- 
lations in analogy with those derived in Appendix [5] f° r 
{ ( n K ) Si } i the average multiplicities of the hadron species 
k emitted by the particular blob i, and the correspond- 
ing second moments, {{n$ n/^) ■ (This approach can 
of course be extended to higher moments.) 

We first note that the canonical partition function for 
N blobs can be expressed recursively in terms of those 
for fewer blobs. For example, with N = N' + (N — N'), 

I S (2N( S ) = J2l s <(2N'Q Is-S'(2(N-N')t s ) . (C4) 

S' 

The inclusive probability for one particular blob to have 
a given strangeness can then be expressed as follows, 

jf-^OSi) = J2 p(Si,...,s n ) 

S2---Sn 

_ / Sl (2&) 7_ Sl (2(JV-l)Ca) 



(C5) 



while the inclusive joint probability for two given blobs 
to have specified strangenesses is of the following form, 



pg- N \s 1> s 2 ) = ]T p(s 1 ,s 2 ,...,s n ) 

S3...SN 

IsMs) IsMs) I- Sl -sMN-2)( s ) 



Io(2N( s ) 



,(C6) 



It is then straightforward to express the average multi- 
plicity of the hadron species n arising from a given blob, 



(l)x(l-JV) 



Ep (i.-*), 

Si 



(S 1 )(n^) Sl • (C7) 



A similar expression holds for the average of any power 
of that multiplicity, (n« ) m . It also readily follows that 
the mixed multiplicity moment for two hadron species 
emitted from the same blob is 

Si 

while the corresponding expression for emission from two 
different blobs is 

\ n K n K> )0 - / , F 12 \ b ll b V\ n K ISA™*' /S 3 ■ 

S1S2 

(C9) 

These expressions allow us to use the canonical partition 
functions for individual blobs obtained in Appendix El 
to calculate the average multiplicities of specific hadron 
species and any desired (co)variances in terms of the cor- 
responding expressions for canonical emission from a sin- 
gle blob. 
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